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Abstract. We introduce a new class of optimization problems called inte- 
ger Minkowski programs. The formulation of such problems involves finitely 
many integer variables and nonlinear constraints involving functionals defined 
on families of discrete or polyhedral sets. We show that, under certain as- 
sumptions, it is possible to reformulate them as integer linear programs, by 
making use of integral generating sets. We then apply this technique to the 
network design problem for fractional and integral flows subject to survivabil- 
ity constraints. 



1. Introduction 

Let us consider the problem of optimal design of a telecommunications network. 
In such a network, digital data are to be sent between pairs of terminal nodes. The 
communicating terminals are, of course, not connected directly to each other but 
only indirectly via a network of switching nodes. In operating and maintaining a 
telecommunications network, several optimization problems need to be solved, each 
on a different time scale. 

For the routing problem, the network is considered fixed, and certain demands 
are given. The direct connections (links) between the switching nodes of a given 
network have certain given transmission capacities, usually measured in units of 
GiBit/s. The objective of the routing problem is to find a way to send the data 
along paths of the network, so that the transmission capacities on the individual 
links are not exceeded. This makes the routing problem a classical problem of 
combinatorial optimization, namely a multicommodity flow problem. However, the 
routing problem usually is complicated by further requirements, like quality-of- 
service (QoS) guarantees. For instance, for voice and multimedia communications, 
it is essential to bound the delays, i.e., the amount of time between sending and 
receiving a data packet, along the paths of a routing for the data of a terminal pair. 
In these situations, it is often also requested to send the data of a terminal pair 
along a unique path of the network, i.e., the flow cannot be split. 

On a larger time scale, of course, the capacities of the links are not fixed. The 
capacities are determined both by the characteristics of the cables and by the 
switching nodes. The cables are replaced very rarely because this is extremely 
expensive; however the capacities of most links can be increased by replacing the 
transmitter-receiver technology installed in the switching nodes. We note that, in 
fiber-optic technology, the installed capacities are usually symmetric. In the capac- 
itated network design problem, or more precisely the capacity expansion problem, a 
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set of (projected) demands are given, which the current network cannot satisfy; the 
objective is to find a set of links and new capacities, such that the demands can be 
satisfied, and such that the installation costs are minimal. On an even larger time 
scale, also the addition of new links, changing the topology of the network, can 
be planned and implemented, to increase the capability of the network. Network 
design problems have been co n sidered in a large body o f literature, see fo r instance 
iBienstock and Gunliikl l|l99fih . iBienstock et~al1 <|l99£l) . lAtamtiirkl l|2002l) and the 
references within. 

In network design problems, one is often concerned with survivability require- 
ments. When one link or several links fail, there should still be a routing that 
satisfies all demands; this requirement is usually called arc survivability. Similarly, 
when a node fails, there should still be a routing that satisfies all demands. How- 
ever, the demands where the failing node is the source or the destination need not be 
considered. This requirement is called node survivability. Network design problems 



for instance iGro 


tschel and Monmal dl99C 


.iGrotschel et al.1 (ll992albJ). IStoerl ( 


l£9j 


Stoer and Dahll 


|l994lll998MLisser et al 


l|l995albh. IBienstock and Muratore 


Noil) 



In the following we will discuss some of the models and solution approaches 
regarding network design problems under survivability constraints. To this end, let 
us fix some notation. Let c denote a vector of installed capacities, and let d denote a 
vector of demands; we shall make the spaces where these vectors live precise in the 
following sections. We shall denote by T c ^d the set of feasible routings. Again, the 
precise space that T c ,d is a subset of will depend on the choice of the formulation; 
w e will discuss it later. 

IStoer and Dahll <ll994h seem to have been the first to address survivability ques- 
tions in capacitated networks; before that paper, research had focused on modeling 
graph-theoretic connectivity constraints. They based their model on the charac- 
terization of capacit ies that a l low a fractional routing of t he demands by so-called 
metric inequalities |lrH . Il97ll lOnaea and Kakushol Il97l() . Therefore their model 
does not contain variables for the no- fault routing and the fault routings; such a 
model is known as a projected formulation. 

There are two shortcomings in this approach. First, allowing fractional routings 
often is too weak a relaxation of the problem. For the existence of an integer 
multicommodity flow, however, there is no known polyhedral description in general. 
Therefore, no projected formulation for the case of integral routings is known in 
the literature. Second, more complicated constraints like those coming from QoS 
guarantees cannot be incorporated into the model of this approach. For instance, 
bounds on the admissible delays along paths of the rou t ing ca nnot be modeled, 
since there are no path variables present. IStoer and Dahll (|l994l) and various other 
authors consider the model of complete rerouting, i.e., in the case of a fault the new 
routing can be completely unrelated to the no-fault routing. 

Other models of survivability require that in the case of a fault only the in- 
terrupted flow needs to be rerouted. There are two main models here. In the 
model of local rerouting, the failure of a link e = (u, v) carrying aggregated flow 
Xg SS in the no-fault routing creates a new demand to route Xg gs units of flow from 
node u to node v in an auxiliary network where the failed link e is not present. 
Local rerouting is also c alled link-based re- establishment in parts of the literature 
()Pi6ro and Medbl I2004T) . In the model of global rerouting, first the flow along 
paths that use the failed link is removed from the no-fault routing; in this way, 
some capacities on links along these paths are released. Then the removed parts 
of the demands are routed through the auxiliary network. Clearly global rerouting 
is more economical than local rerouting, since the released capacities can be used 
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for the fault routings . G lobal rerouting is also called path-based re- establishment 
l|Pioro and Medhil . I2004L 

Contribution of this paper. In this paper we investigate the structure of the 
family of feasible-routing sets T c .d when the installed capacities c and the required 
demands d vary. In our mathematical model, both the capacities and the demands 
can take arbitrary (non-negative) integer values. Thus the family of feasible-routing 
sets is infinite. The key observation in order to be able to deal with this infinite 
family is that Minkowski sums of these sets play an important role. We will show 
that in general the feasible-routing sets satisfy 

^cKd 1 + F c 2_ d 2 C T c i +c 2_ d i +d i\ (1) 

then we consider the pairs (c 1 ,^ 1 ) and (c 2 ,d 2 ) where the inclusion is in fact an 
equality. 

First, it is possible to address the question whether a capacity installation c is 
feasible for a demand vector d. Whenever there is a pair (c 1 ,^ 1 ) and (c 2 ,d 2 ) such 
that c = c 1 + c 2 and d = d 1 + d 2 and in £[J equality holds, then c is feasible for d 
if and only if c 1 is feasible for d 1 and c 2 is feasible for d 2 . 

Second, it is possible to address survivability questions as well. As a first exam- 
ple, let a be an arbitrary link (arc) of the network. We will define a real-valued 
functional g a on the family of feasible-routing sets by defining g a {J- c .d) to be the 
minimum amount of flow carried on a in any routing of J- C ,d- Then the network with 
capacity installation c and demand d is arc-survivable (in the model of complete 
rerouting) if and only if g a (J- C: d) = for all arcs a. It turns out the functional g a is 
Minkowski- additive. This implies that, if equality holds in |T]). then (c, d) defines a 
survivable network if and only if (c 1 , d 1 ) and (c 2 , d 2 ) do. 

Outline. The outline of the paper is as follows. In lsection 21 we introduce a class 
of mathematical optimization problems called ( mixed- )integer Minkowski programs 
as our general model; their formulation includes constraints involving non-linear 
functionals like g a . We prove a general reformulation theorem that asserts that 
mixed- integer Minkowski programs can be reformulated as integer linear programs, 
whenever there exists a suitable finite generating set. 

In Isection 3l and Iscction 4l we show that in two important settings, related to 
fractional and integral routing, there exists a finite generating set for the families 
of feasible- routing sets T c ,d- In lsection 51 we make the formulations of the network 
design problems precise and discuss the structure of the generating sets that arise. 
This enables us to reformulate the network design problem as an integer linear 
program. 

In lsection 61 wc address questions of survivable network design in the model of 
complete rerouting both for fractional and integral flows. Again the survivability 
conditions can be modeled using Minkowski-additive functionals in the framework 
of integer Minkowski programs. The general reformulation theorem, applied to 
this formulation, then gives a new formulation of the survivable network design 
problem as an integer linear program. The resulting formulation as an integer 
linear program does not contain variables for the no-fault routing and the routings 
in the individual failure scenarios. Hence, the formulation can be regarded as a 
projected formulation. This is particularly relevant in the case of integral routings, 
where no projected formulation was previously known. 

A prospective advantage of the proposed formulation of the survivable network 
design problem is that the formulation only grows moderately with the number of 
failure scenarios. Indeed, only one linear inequality and no extra variable is needed 
for each failure scenario. The number of variables only depends on the topology 
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of the network. Moreover, the reformulation only needs to be computed once for 
a given network topology; it can then be used for solving the survivable network 
design problem for arbitrary demands. On the other hand, the usual non-projected 
formulations contain one set of routing variables for every failure scenario. 

2. Integer Minkowski programs 

2.1. Definitions and first results. In this section we define (mixed-)intcgcr 
Minkowski programs and show that under certain conditions there exist reformu- 
lations as integer linear programs. (Mixcd-)integer Minkowski programs consist of 
integer variables, a linear objective function and inequality constraints involving a 
certain class of nonlinear functions. We first remind the reader of the definition of 
Minkowski sums. 
Definition 2.1. 

(i) Let Ax,A 2 C R d be sets. The Minkowski sum (pointwise sum) of A\ and A 2 
is defined as 

Ai + A 2 = { x 1 + x 2 : x 1 G Ax, x 2 G A 2 }. 

(ii) Let A C R d be a set and let k G Z + be a non-negative integer. We will denote 
by kA the set obtained by taking the Minkowski sum of k copies of the set A, 

kA= {x 1 + --- + x k : x 1 ,...,x k G A}. (2) 

In the same way, we will use the notation Yliel kiAi. 

Note that when one of the summands is the empty set, then the Minkowski sum 
is empty too. 

Remark 2.2. Frequently, Minkowski sums are considered in the context of convex 
sets only. For convex sets, the notation kA is usually defined for arbitrary scalars 
k eR as 

{ kx : x G A } ; 

for k G Z + this gives the same set as in lf2"jl. However, when A is not a convex set, 
the set kA as defined in J2J is not necessarily the same as the set { kx : x G A } . 

We will consider a class of nonlinear functions that are obtained in a two-stage 
process. In the first stage, we define a set-valued mapping A: Z™ — > 2 R + that is 
superadditive: 

Definition 2.3. A set-valued mapping A:~L n — > 2 R + is called superadditive if it 
satisfies 

A(z x + z 2 ) D Aiz 1 ) + A(z 2 ) for Zl ,z 2 eZ n , (3) 
where + in the right-hand side denotes a Minkowski sum. 

In the second stage, we define a functional k on the range of A with the following 
properties: 

Definition 2.4. Let n: 2 K + D A(Z n ) -> R = R U {+00} be a functional with 
«(0) = +00 and A,B,C G A{Z n ). 

(i) k is called Minkowski-additive if C = A + B implies k(C) = k(A) + n(B). 

(ii) k is called antitone if A C B implies k(B) < k(A). 

With these preparations, we are now in the position to define (mixed-)integer 
Minkowski programs. 
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Definition 2.5. A (mixed-)integer Minkowski program is an optimization problem 
of the following structure: 

min W T z 
s.t. A{z)^$ 

B ■ z = v (4) 
k 1 {A{z)) + hjz < Si fori=l,...,k 

z e T\ 

where W G R n is a cost vector, A: Z" — > 2 K + is a superadditive set-valued mapping, 
B G Z nxn a matrix, v G Z™, n l : A(Z n ) — > R /or i = 1, . . . , k is a Minkowski- additive 
and antitone functional, and hi G R™ /or i = 1, . . . , k. 

In the mixed-integer case, we have A(z) C for all z G Z™ , whereas in the 
integer case A(z) C Z^ /or aZZ z G Z n . 

The following theorem will concretize the conditions under which we can refor- 
mulate a (mixed-)integer Minkowski program as an integer linear program. 

Theorem 2.6. Let {z^}j C Z™ with A(z^) ^ for all j € J be a finite generating 
set of the family of sets A{z), i.e., for z € Z™ there is a representation 

z = Xj z J where Xj G Z + , 
je>7 

ITien the ( mixed- Jinteger Minkowski program (@J) can be reformulated as an integer 
linear program: 

min ^2(W T z j )Xj 

s.t. E (Bz j )Xj = v 

jej (5) 

E [n i (A(z i )) + hjz 1 } Xj <5 { fori = l,...,k 

Xj G Z+ for j G J 

Proof. Let A*, j e J be a feasible solution of the integer program |J5J|. We need 
to prove that z* = J2jej & * s a feasible solution of the corresponding (mixed-) 
integer Minkowski program (0}. Relation © yields: 

A(z*) =a(^2 A i * j ) 2 E A i ( 6 ) 

and as A{P) ^ this implies A(z*) ^ 0. We also have B ■ z* = v. It remains 
to prove that the linear constraints of the (mixed-) integer Minkowski program are 
satisfied. Let i G {1, . . . , k}. As k 1 is an antitone functional for all i = 1, . . . , k and 
as © holds, we have: 

k 1 (A(z*)) + hfz* < E A* n l {A{z J )) + E Xj hfz j < Si. 

This proves the feasibility of z* = K & f° r 0- 

jeJ 
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On the other hand, let z be a feasible solution of the (mixed-)integer Minkowski 
program (@J. We know that there is a decomposition z = J2jej Pj ^ such, that 

where f3j G Z + for all j G J\ We shall prove that j e J, is a feasible solution 
of the integer linear program JSJ. Surely we have B ■ z = 5Zj gi y (Bz^)f3j = v. Let 
i G {1, . . . , k} be fixed. As k j is a Minkowski-additive functional for alH = 1, . . . , k 
we have 

Si > k 1 (A(z)) + h[z=Yl & l Ki ( A (z j )) + hjz j }. 

Therefore. f3j, j G J is a feasible solution of the integer linear program JHJ). 

Clearly also the objective value of the integer linear program {JjJ is equal to the 
objective value of the corresponding (mixed-)integer Minkowski program Q for 
every cost vector W G R n . □ 

2.2. Families of truncated sets and truncated integer Minkowski pro- 
grams. The results of the previous subsection may be improved in the sense that 
it is possible to formulate integer Minkowski programs on families of truncated sets 
as integer linear program. This result is important for our application in network 
design because there are usually lots of circulations of flow as feasible solutions 
for the multicommodity flow problem which are not used for a "regular" objective 
function and therefore should be cut away. We give the definition of a truncated 
set first and continue giving the results for truncated (mixed-)integer Minkowski 
programs. 

Definition 2.7. Let u,v G R n be vectors. We say that u reduces v and denote 
uC« if u^v^ > and \u^ \ < \v^'\ for all components i = 1, . . . ,n. 

Definition 2.8. Let C C R n (or C C Z n ) be a set of vectors. Let A: Z" -> 2 R + be 

a set-valued mapping which is super- additive and z G Z™. We call 

tv c {A{z)) = {x G A(z) :$c,eC with c C x} 
the truncated set corresponding to z and A with respect to C. 

The definition of truncated sets is illustrated in Figure ^ f° r a discrete set of 
vectors and a polytope. Both are truncated by a fixed set C of two vectors. 

Definition 2.9. A truncated set trc(^4(z)) is called indecomposable if there is no 
^ z 1 ,z 2 e Z™ with 

tr c (A{z)) C tT C (A{z 1 ))+tr c (A(z 2 )). 

Lemma 2.10. Lf trc(A(z)) is indecomposable then A{z) is indecomposable, too. 

Proof. We have tic(A(z)) C A(z). Suppose there is a decomposition A(z) = 
A(zi)+A(z2). Let x G tvc(A(z)) C A(z), then there is a decomposition x = x\ +x 2 
with Xi G A(zi), i = 1,2. In particular: a;,; C a; (because Xj G Suppose there 

is c G C with (w.l.o.g.) c C zi, implying c □ 2, which is a contradiction. Thus 
2j G trc(^4(zi)) anc i we have tic(A(z)) C trc(^4(zi)) +trc(^4(z2))j contradicting the 
indccomposability of trc(^l(z)). □ 

The previous lemma shows that there is a finite generating system of families 
of truncated sets provided there is such a finite generating system for the corre- 
sponding family of sets. We go on showing that under certain conditions it is 
possible to reformulate integer Minkowski programs on truncated sets as integer 
linear programs. 
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Figure 1 . Truncation for a discrete set of vectors and a polytope 

Furthermore, let W G R" denote a vector of costs, B G j nx ™ a matrix, »eZ", 
k 1 : A(Z n ) — ► R for i = 1, . . . , k is a Minkowski-additive and antitone functional, 
and /ij £ R" for i = 1, . . . , k. We consider the following truncated integer Minkowski 
program: 

min W T z 
s.t. tr c (A(z)) ^ 

B ■ z ~ v (7) 
K*(tr c (.4(z))) + /if 2 < 5j for i = 1, . . . , k 

z e z". 

Lemma 2.11. Let {z~j}jej denote the finite generating system of a family of trun- 
cated sets, i.e., for z G Z n there are Xj G Z+, 7 G J7\ such that z — ^2 ^jZj 

and 

tx c {A{z)) cJ^Aj-treC^)) 
ie.7 

7/K I (trc(^(z))) = /or all z G Z n and i = 1, . . . , k then we can reformulate 

the truncated integer Minkowski program J7J as integer linear program: 

min ^{W T z r )X ] 

s.t. ]T (B^)Ai = » 

je.7 (8) 

53 [k 4 (tr c ) ) ) + hf z*] Xj<S t fori = l,...,k 

Xj G Z+ for j ej 

Proof. Let z be a feasible solution of the truncated integer Minkowski program (JZJ . 
Then there are (3j G Z + such that z = ftj % 3 an d 
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It is clear that ^ (Bz^Xj = v holds and as the functions n % are antitone for 
i = 1, . . . , k we have: 

St > n\tT C (A(z))) + hfz>Y, [K l (tr c (A(zi))) + hfz j ]pj 

for i = 1, . . . , k. This means that f3j £ Z+ for j £ is a feasible solution of the 
integer linear program (JHJ). 

It remains to show the other direction. Let Xj e Z + for j e J be a feasible 
solution of the integer linear program (JHJ - We set z := Aj f J ' and show that z 

is a feasible solution of the truncated integer Minkowski program (JJJ. We have: 

v = Y i iP#)\ j = B ■ z. 

On the other hand, as z = Aj 2 J we have ^4(z) 2 X) Aj ^(z- 7 ). This implies: 

jej jeJ 

5 i >Y d [^{^c{A(z i ))) + hJz j ]\ j 

= J2[>j(Az j )) + hfz3]X j 
jej 

> K l {A{z))+hf Z 

= K i (tv e (A(z))) + hfz, 

for i = 1, . . . , k. This means that z is a feasible solution of the truncated integer 
Minkowski program Q. □ 

In the following sections, we will investigate some finitely generated families of 
sets A{z) and the structures of the associated truncated (mixed-) integer Minkowski 
programs. 

3. The integer case: Atomic fibers 

In this section we will investigate the situation when the mapping A generates a 
family of discrete sets. In our application, the discrete sets will be sets of feasible 
integral routings. To be more precise, let A E x dxn be a matrix, and let the 
mapping A: Z d — > 2 R + be defined by 

A(b) = P I A . h = {z : Az = b, z E T\ }. (9) 

T he set P\ h is known as the f iber of b under the linear map Ja ■ Z™ —> Z d , x i— > Ax. 

lEisenschmidt et all l)2006a() considered the family of fibers of a fixed matrix A £ 
Z dxn , when the right-hand side vector b £ Z d varies. They established a theory of 
Minkowski decomposition for fibers; in this context, the non-decomposable fibers 
are called atomic: 

Definition 3.1. We call a fiber P\ h atomic, ifP^b ~ ^Ab +^A& implies b = b\ 
or b = b-2- 

For our purposes, we will need to consider a slight generalization of the notion of 
atomic fibers, where we restrict the set of right-hand sides from ~L d to some subset. 
As a motivating example, consider flow-conservation constraints that appear in 
standard node- arc formulations. These constraints are linear equations A^.z = bi 
with a right-hand side bi — 0. It is desirable that the flow-conservation constraints 
also hold in the Minkowski summands P A bl , b2 of a set Pj± ^ of feasible integral 
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routings. Therefore, we wish to restrict the decomposition into atomic fibers such 
that bi = for all components corresponding to the flow-conservation constraints. 

This example suggests to restrict the right-hand sides to a sublattice of Z d , but 
later we will see that it is also useful to consider restriction to a submonoid of Z d . 

Definition 3.2. Let (M, +) with M C Z d be a monoid, (A, +) a lattice with A C 
Z d , and let A £ Z dxn be a matrix. 

(1) A fiber P A b with b £ A is called atomic w.r.t. A if there is no decomposition 
b = b 1 +b 2 with P\ b = Pj bl + P' A b2 , where b 1 ,b 2 £ A. 

(2) A fiber P A b with b £ M is called atomic w.r.t. M if there is no decompo- 
sition b = b 1 +b 2 with P ! Ab = P ! AM + P A b 2, where b 1 , b 2 £ M. 

Once again we ref er to the theory of indecomp osable fibers w.r.t. monoids and 
lattices developed by lEisenschmidt et "all l|2006al) and recall a result which is im- 
portant for our analysis in this paper: 

Lemma 3.3. Let M C Z" be a finitely generated monoid and let A £ Z dxn be a 
matrix. There are only finitely many fibers of the matrix A that are atomic w.r.t. M. 

Note that for M = (±A±, . . . , ±A n ), where Aj, i = 1, . . . ,n, denote the columns 
of matrix A, the above lemma shows that there are only finitely many atomic fibers 
of a matrix A. 

3.1. Truncated fibers. The usual formulations of the routing problem (multi- 
commodity flow problem) include routings that contain certain flow circulations as 
feasible solutions. When a useful ("regular") objective function is chosen, an opti- 
mal solution to the routing problem will never contain a flow circulation, when we 
consider the routing commodity by commodity. The reason is that a routing with 
circulations will always be dominated by the routing where the circulations have 
been removed. For our purposes, however, it makes sense to cut away circulations 
explicitly; we will see later that it can significantly simplify the computations. 
This application gives rise to the notion of truncated fibers. 

Definition 3.4. Let C = {ci, . . . , c s } with Ci £ ZI? be a finite set of vectors. Let 
P Ab be a fiber of a matrix A £ Z dxn . We call 

tr c(PA.b) = {z£ P{,b ■■ $c£C with c C z} 
the truncated fiber of A with right-hand side b with respect to C. 

Definition 3.5. A truncated fiber tie (P A b ) is called indecomposable (w.r.t. a monoid 
M ) if there is no b\ , 62 7^ (b\ , 62 £ M) with b = b\ + 62 and 

tr c (Pj >6 ) Ctr c (Pj >6i ) +tr c (Pl i6 J. 

Corollary 3.6. Ifttc(PA,b) * s indecomposable (w.r.t. a monoid M), then P Ab is 
atomic (w.r.t. the monoid), too. 

Proof. This is a direct consequence of Lemma 12.101 □ 

This leads to the main result of this section: We can reformulate the integer 
Minkowski program on truncated fibers as an integer linear program. Consider the 
integer Minkowski program: 

min W T b 

s.t. tr c (Pl i6 )^0 

B-b = u (10) 

^(trciPi^) + hfb<S l for i = 1, . . . , k 

b £ Z d , 
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where W £ R d is a vector of costs, B £ / dxd j u g z d , k 1 antitone and Minkowski- 
additive functionals, hi £ R d and Si £ R for i = 1, . . . , k. Then we have: 

Corollary 3.7. If ^(ti c [Pi jb ) = for all b £ T d and i = 1, . . . ,k then we 

can reformulate the integer Minkowski program (jlOfl as an integer linear program: 

min ^(W T V)\ 3 

s.t. >T (ByjAj- = u 

jej (11) 

]T [«*(tr c (Fi >5J ) + fcf &>'] Aj < <5i /or i = 1, . . . , k 

Xj £ Z+ for j £ J 

Proof. This is a direct consequence of Lemma 12.111 □ 

3.2. Projections of fibers. The two basic formulations for the routing problem 
are the node-arc formulation and the path formulation. It turns out that the 
node-arc formulation essentially is a projection of the path formulation. By the 
term projection we mean linear, integral transformations of problems from a high 
dimensional space of variables to a space of variables of lower dimension. 

Because of this relation between the node-arc and path formulations, we also 
need to consider projections of (atomic) fibers. We consider the fibers of the inte- 
gral matrix A £ Z dxn under a projection described by a matrix II 6 Z™ x ". We 
remark that projections of fibers are not nec essarily fibers of some matrix A; they 
can have a more complicated structure, see IWilliamsl l|l992j) . In our application, 
however, the following setting is general enough. We suppose that there is an inte- 
gral transformation matrix N £ 1° xd for the set of right-hand sides b £ Z d , with 
the columns of N being linearly independent. Furthermore we suppose that there is 
an an integral matrix A £ Z cxm such that TV^! = AH. Then we have the following 
property: 

II • P\ h = II • {z £ ~l\ : Az = b} = II • {z e Z™ : NAz = Nb} 

= n-{zeZ" + : AYlz = Nb} (12) 
= {y £ Z™ :Ay = b} = P^, 
where b = Nb. Equation (|12|) gives the following lemma. 

Lemma 3.8. If P\ b is decomposable thenH-P^ b = P^ g is decomposable w.r.t. the 
monoid B = {b = N ■ b : b £ AT%}. 

Proof. This is a direct consequence of the linearity the projection LI: 

pi tl = iL.pi !b = n-(p{ M +Pi ib2 ) 

= LI • P{ M + LI • Pi M 

_ pi _ i pi _ 

□ 

An analogous assertion is true for truncated fibers. 

Lemma 3.9. Let C = {ci, . . . ,ct} C Z™ a finite set of vectors and letT> C Z™ such 
that A ■ d £ B for all d £ T>, where B is the monoid from Lemma \3.lA For all c £ C 
let there exist d £ D with d □ 7r(c). Then we have: If trc(Pj^ b ) is decomposable, 
then tvT>(Pj^ b ) is decomposable w.r.t. B. 
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Proof. As a first step we will show that tr c (P| £ ) = tr p (II • P^ b ) C II • tr c (P{ b ). 
To this aim let x G trx>(II ■ P^ b ). Then there is y G P\ b with II ■ y = x. Suppose 
there is c G C with c C y. Then there is c? G V with dCII-cCII-y = a; which 
contradicts the fact that x G trx>(P^ b ). Thus y G trc(P^ b ) and consequently 
a; = II • y G II • trc(P^ b ). This proves our first claim. Now let trc(P^ b ) be 
decomposable, i.e., there are vectors bi, hi with tre(P^ b ) Q trc(P^ b ) + trc(P4 b )■ 
By linearity of projections we have 

n ■ tr c (p| >6 ) c n • tr c (p{ m ) + n • tr c (Pl ib2 ). 

We will show that trz»(Pj b ) C trz,(Pj 5 ) + tr^PjgJ. Let x G tr c (Pj b ) C 
tt-tr c (P{ b )). Thus there is a representation x = x\ + X2 with xi G LT • trc (PjJ b ). 
We will show that Xi G trx>(P| b .) for i = 1,2. Suppose not and let w.l.o.g. x\ ^ 
trp(P^ bi ). Then there is d G I? with dCxi. But then g? C xi + X2 = x because 
Xi G Z™. This contradicts the fact that x G trx>(P^ b ) and concludes the proof. □ 

This implies: If we have a finite generating system for the fibers P\ b (or 
tre(P_4 b )) then we have a finite generating system for the projected fibers P\ b 
(or tru(P^ b )) w.r.t. the monoid B. Therefore we may reformulate the integer 
Minkowski program on the projected fibers as an integer linear program, because 
we may do so for the original fibers. 

4. The mixed-integer case: Indecomposable Polytopes 

4.1. The general case. In this section we will consider finitely generated fami- 
lies of polyhedral s ets depending on integral vectors. We refer to the notation of 
iHenk et alJ ^OO.^ . IKoddI l|20021) . Let W G Z mx " be a fixed but arbitrary integral 
matrix with row vectors w l G Z™, 1 < i < m. We assume that 

pos{w\...,w m } = R n , 

where pos denotes the positive hull. Thus for every u G R m the set P u = {y G R™ : 
Wy < u} is a polytope. We are interested in the set of all nonempty polytopes 
with integral right-hand side arising in this way. We set: 

U{W) = {u G R m : P u ± 0}, 

and consider U{W) n Z m . 

Definition 4.1. A polytope P z , z G U(W)C\~L m is called integrally decomposable if 
there exist P z i , P z 2 not homothetic to P z such that P z = P z i + P z 2 and z = z 1 + z 2 , 
z l G U{W) n Z m . P z is called integrally decomposable otherwise. 

We have the following result in IHenk et alJ l|2003l) : 

Theorem 4.2. There exist finitely many vectors h x ,...,h k G U(W) n Z™ such 
that for every polytope P z , z G U{W) R Z m 7 there exist h n , . . . , h j2m - 2 ~ n and non- 
negative integers Xj 1 , . . . , A 3 - 2m _ 2 _ n such that 

2m — 2 — n 2m — 2 — n 

Pz = ^jiPfiii an d z= ^jiti**- 

i=l i=l 

This means in particular, that it is possible to model a mixed-integer Minkowski 
program on polytopes with integral right-hand side as an integer linear program 
according to Theorem 12. Gl 
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4.2. Restricted right-hand sides. In our application to network design the poly- 
topes P z introduced in the previous subsection will represent the sets of feasible 
vectors of flow of particular networks. Therefore it makes sense to claim non- 
negativity of the points in a polytope P z . We assume that our matrix W is of a 
special structure: W = (W, — Id) T G Z' m+n ' x ™ and we restrict our attention to the 
following lattice of right-hand sides: 

* = {u = («, 0) T G Z m+n : u G Z m }. 

For u6f we have: 

Pu = {y G R n : Wy < u} = {y G R™ : Wy < «}. 

As in the integral case of atomic fibers we are interested in sublattices of iff, be- 
cause of flow-conservation constraints or non-negativity constraints for demands 
and capacities for example. 

Definition 4.3. Let A C $ denote a lattice and M C iff denote a monoid. 

(i) A polytope P u with u G A is called integrally indecomposable w.r.t. A if there 
are no vectors ^ u\,ui 6 A with P u = P Ul + P U2 and u = u\ + u%. 

(ii) A polytope P u with u G M is called integrally indecomposable w.r.t. M if 
there are no vectors ^ U\, U2 G M with P u = P Ul + P U2 and u = Ui + it2- 

It was already proved in iHenk et alJ l)2003|) . Remark 3.1, that there are only 
finitely many integrally indecomposable polytopes w.r.t. a sublattice A of Z m . We 
may extend this result to monoids under certain conditions. Let M be a monoid 
which is finitely generated, i.e., which is generated by mi, . . . , m t £ M and let Am 
the lattice generated by mi , . . . , m t . 

Lemma 4.4. Let M be a monoid and let W~L\ n Am Q M . Then there are only 
finitely many integrally indecomposable polytopes w.r.t. M . 

Proof. We know that there are only finitely many integrally indecomposable poly- 
topes P z w.r.t. Am- As M C Am wc know, that if z G M and P z integrally 
indecomposable w.r.t. Am then P z is integrally indecomposable w.r.t. M. On the 
other hand, let P z = P Zl +P Z2 with z = Z1+Z2, z, Z\, Z2 G Am and P z ,P Zi ^ 0. Then 
z, Zi G WZ n , i = 1,2. As WZ n n Am = WT% n A M Q M, we have z,z u z 2 G M. 
This means that the polytopes which arc integrally indecomposable w.r.t. Am are 
integrally indecomposable w.r.t. M and vice-versa. This implies the finiteness of 
the number of integrally indecomposable polytopes w.r.t. M. □ 

Corollary 4.5. Let C be a set of vectors (not necessarily finite). Consider the 
truncated polytopes ^ trc(P«) = {y G R" : Wy < u, pc G C with c\—y}. There is 
a finite generating set for these truncated polytopes. 

Proof. This is a direct consequence of Lemma I2.1UI □ 

Corollary 4.6. One may reformulate integer Minkowski programs on (truncated) 
polytopes with integral right-hand side as integer linear programs. 

Proof. This comes from Theorem 12 . 1 II and Thcorcm l2.6l □ 

5. Network design problems and atomic fibers 

In this section, we will treat the network design problem in terms of the analysis 
of atomic fibers. The outline of this section is as follows: The first two subsec- 
tions will introduce two formulations of the network design problem, the node-arc 
and the path-cycle formulation. The following subsection will give the connection 
between the path-cycle and the node-arc formulation and introduce the notion of ir- 
reducible networks. The following two subsections will give the connection between 
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the irreducible networks and the atomic fibers for both our formulations. The last 
subsection will finally present the network design problem as an integer Minkowski 
program and its reformulations as an integer linear program. 

The reason for studying both formulations is the following. The path formulation 
is easier to study than the node-arc formulation: In the path formulation, ordinary 
atomic fibers appear, whereas in the node-arc formulation, we need to consider 
atomic fibers with respect to a monoid of feasible right-hand side vectors. On the 
other hand, the node-arc formulation is essentially a projection of the path-cycle 
formulation. Therefore, the set of atomic fibers is in general much smaller in the 
node-arc formulation. Consequently, the computation of the atomic fibers is more 
efficient, and also the reformulation as an integer linear program is more compact 
when we start with the node-arc formulation. 

Let G = (V, A) be the supply digraph, which is connected. Each arc a G A of 

\A\ 

the graph has a capacity c a G Z + such that c G ~Lj_ and we are given a demand 
vector d = (di, . . . , dk), di G Z + Vi G {1, . . . , k}. For each commodity I = 1, . . . , k 
let si G V be its source, ti G V its sink and II; the set of all paths from si to ti . Such 
a setting is called a multicommodity network. We will denote it by N = (V, A, d, c) . 

5.1. Node-arc formulation. In this subsection, we will give an exact definition 
of the problems we consider in this paper. Let /• G Z+ be the part of the flow 
of commodity I which uses arc a. The problem of finding a feasible flow w.r.t. the 
capacity vector c and the demand vector d can be formulated as follows: 

k 

J2f a +*a =c a Vac A (13a) 

i=i 

E E fa = d l VZ = l,...,fc (13b) 

aeS+(s,) aeS-(st) 

E fa' E fa =0 VsenKt,}, Vi = l,...,A (13c) 

aeS+(x) aeS-(x) 

f a ,s a eZ+ Va€A, VZ = l,...,fc, (13d) 

where s a , a G A, denote the slack variables. We denote by C the matrix corre- 
sponding to the system of equations Ijl3|l . 

Definition 5.1. Let G = (V,A) be a digraph, d = (di,...,dk) G J-\ a given 
demand vectors and c G ~L + a given capacity vector. Let K a G R+ denote the costs 
of routing one unit of commodity I through arc a. The program 

min K T f 

s.t. C-(f,s) T = (c,d,0) T (14) 

f a ,s a ez + VoeA,VJ = i,...,* 

is called the multicommodity flow problem for integer flows with respect to the 
node- arc formulation. 

In the multicommodity flow problem, we are given a demand vector and a capacity 
vector. We want to find a routing of the demands, which is minimal w.r.t. the costs 
of routing flow. In the network design problem, we are given a vector of demand 
and want to find a capacity vector minimizing installation costs, such that we can 
find a routing for the demands. 

Definition 5.2. Let G = (V,A) be a digraph, d = (di,...,dk) G Z5_ a given 
demand vector. Let additionally W a , a G A, denote the costs of installing one unit 
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of capacity on arc a. The program 
min W T c 

s.t. C-(f,s) T = (c,d,0) T 

fls a eT+ VoGA,V/ = l,...,A 1 j 

c a g Z + Va e ^4 

?s called the network design problem for integer flows with respect to the node- arc 
formulation. 

5.2. Path-cycle formulation. Let W be the set of directed cycles in G = (V,A). 
y l p denotes the part of commodity I routed on path p, p G 11; . y l w denotes the part 
of commodity /, which circulates on cycle w G W. The problem of finding a feasible 
flow w.r.t. the capacity vector c and the demand vector d can be formulated as 
follows: 

E ( E + E f«) + s « = c « ^ aeA ( 16a ) 

J2v l P =di VIG{1,...,*} (16b) 

pen, 

4,»Lez + Vpen ;i v/e{i,...,fc}, (16c) 

where s , agi denote the slack variables. 

As in the previous subsection we define the network design problem and the mul- 
ticommodity flow problem for integer flows although this time w.r.t. the path-cycle 
formulation. With the notation of the previous subsection, we have the following 
formulation of the multicommodity flow problem: 

mm £ e Rl a( E y l P + E vi, 

aaAi=i \pen, wew / (17) 

a£p a£w 

s.t. est, 

and the following formulation of the network design problem: 

min Vy T c 

s.t. eg (is) 

c a G Z + Va e A. 

Now we look at the multicommodity flow problem: As we have non-negative costs 
associated with the arcs (K G R^' fc ), there will be no cycle flow in the optimal 
solution, i.e., y[ = Vw G W, VI G {1, . . . , k}. So we can eliminate the cycle- flow 
variables y w . For the network design problem, too, we can eliminate the cycle- flow 
variables y l w . We obtain the following formulation: 

E^ = d i VZe{l,...,fc} (19a) 
pen, 

k 

E E y l p +Sa ^ Ca Va e A ( 19b ) 

i=i pen, 

a£p 



2/;gZ+ VpeH.VJ €{!,...,&} (19c) 



4 

s a G Z+ Va G A, (19d) 

where s a , a G A, denote the slack variables. 
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Let D denote the matrix corresponding to the left-hand side of formulation (|19|l . 
We obtain a new formulation of the multicommodity flow problem 

mm E E K( E v l P ) 

aeAl=l p£U t 

s.t. D-(y,s) T = (d,c) T ( 2 °) 

y l p eZ+ ypGU h yi = l,...,k 

s a 6Z + Va e A 

and a new formulation of the network design problem 

min W T c 

s.t. D-(y,s) T = (d,c) 

y l p eZ+ Vpen i VZ = l,...,fc 

c a , s a e Z + Va e A. 



(21) 



This formulation is called the pa£/i formulation. 

5.3. The node-arc formulation as a projection of the path-cycle formula- 
tion. It is clear that one can convert feasible solutions of the multicommodity flow 
problem in the path-cycle formulation to to feasible solutions of the node-arc formu- 
lation via a projection. Let y be a feasible integer solution of the multicommodity 
flow problem with respect to the path-cycle formulation. We set: 

t=E^+E^ Va e A, Vie {l,. ..,*}. 

peri; wew 

Then, / is a feasible integer solution of the multicommodity flow problem with 
respect to the node-arc formulation. Indeed, we have integrality of the components 
of / because y is integer. The capacity constraints are respected because: 

E^ + So = E(E^ + E y l w)+ s ° = c a VOGA 

1=1 1=1 S>erii wew ' 

It remains to check whether the flow-conservation constraints are respected. To 
this aim, we look at a path p = {si,xi, . . . ,x n , ti) with y l p > 0. The flow on path 
p respects the flow-conservation constraints for every x% G p. This observation is 
true for every path p £ II. The same assertion is valid for all cycles w with y l w > 0. 
Now we look at some node Xi G V\{si,ti}. We have: 

E fa- E f a = ° VZe{l,...,fc}. 

a£S+(xi) a£S-(xi) 

Analogous arguments yield: 

E fa' E fa = d l VZG{l,...,fc}. 

aeS+(s t ) aes-(s,) 

Therefore the path-cycle flow y determines the node-arc flow / uniquely. 

Also, one can convert feasible solutions of the multicommodity flow problem in 
node-arc formulation to feasible solutions in path-cycle formulation. 

Lemma 5.3 fsee lAhuia et all l|l993|) . Theorem 3.5). Every non-negative arc flow 
f can be represented as a path and cycle flow y (though not necessarily uniquely). 

However, the solution is not uniquely determined in general. 
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1 2 




3 4 

Figure 2 . A network with non-unique path-flow 

Example 5.4. To see an example of this non-uniqueness consider the one-commodity- 
digraph in Figure |3 There are 4 paths from the source s to the sink t: p\ = {1,2}, 
P2 — {1, 4}, = {3, 2} and p4 = {3, 4}. We consider the arc-flow vector / = (1, 1, 1, 1), 
i.e., one unit of flow on all arcs. This flow can be represented in two ways as path-flow: 

y 1 = (yp^vL'Vps^vD = (ijO.o. *) or y 2 = {vlx,vl 2 ,yl 3 ^vli) = (°> *> l M- 

5.4. Irreducible networks in the node-arc and path formulation. We will 
now introduce the notion of irreducibility of networks. To this aim, we consider the 
set of feasible solutions of the multicommodity flow problem in the node-arc and in 
the path formulation. 

Definition 5.5. Let G = (V, A) be a digraph, d £ Z+ the demands and c € z!^' the 
capacities. We denote by T v the set of feasible integer solutions of the multicom- 
modity flow problem on N = (V, A, d, c) w.r.t. the path formulation. Analogously, 
we denote by JT NA the set of feasible integer solutions w.r.t. the node-arc formula- 
tion. 

Definition 5.6. Let G = (V, A) be a digraph, d £ Z^. the demands and c € zi_ ' the 

capacities. We denote by TN V the set of all non-cyclic feasible integer solutions 
of the multicommodity flow problem on N = (V, A, d, c) . This means: If y G J-N F , 
then y contains no circulation of flow. J 7 N NA is defined analogously. Of course we 
have the following relations: 

TN P C T p TN NA C ^ NA . 

To illustrate these definitions, we consider the following example. 
Example 5.7. Let G = (V,A) the digraph in Figure^l 




3 

Figure 3. A digraph with possible circulations 

We have 4 paths in this example: pi — {1,4}, p2 = {1,2,6}, ps = {5,6} and p± = 
{5,3,4}. Let d = 2 and c = (1, 1,1, 1, 1, 1). Then y = {y P1 ,yp 2 ,yp 3 ,y P4 ) = (0, 1,0, 1) is a 
feasible integer solution of the multicommodity flow problem w.r.t. the path formulation. 
This means y 6 T p . But y contains a circulation of flow on arcs 2 and 3. Therefore it is 
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not contained in TN P . The arc flow / which is determined by y is contained in J rNA but 
not in FN NA . 

Now we can give the definition of irreducibility of networks. 

Definition 5.8. Let G = (V, A) be a digraph, d G ~l\ the demand and c £ the 
arc-capacity of a network. A decomposition w.r.t. the path formulation (w.r.t. the 
node-arc formulation) of the network N = (V, A, d, c) is given by a decomposition 
of the capacity vector and the demand vector c = c\ + c%, C\,Ci £ d = 

di + c?2, G?i,d2 £ Z5_, such that the set of non-cyclic feasible integer solutions of 
the multicommodity flow problem can be obtained as the Minkowski- sum: 
FNl + TNl (TN na = TN^ A + TN^ A ), where FN* (TN^ A ) is the set of 
non-cyclic feasible integer solutions of the multicommodity flow problem on Ni — 
{V, A, di, Ci), i = 1,2. // the network can not be decomposed in this way, it is 
called irreducible or indecomposable w.r.t. the path formulation (w.r.t. the node- 
arc formulation). 

The notion of irreducibility depends on the formulation of the multicommodity 
flow problem, as the following example illustrates. 

Example 5.9. Let us consider the network in Figure |1J with all the arc capacities equal 
to 1. 



1 1 




1 1 



Figure 4. A network with demand d = 2 

We have already considered this example in Figure |21 The set of solutions of the 
multicommodity flow problem w.r.t. the node-arc formulation is a singleton: ^ NA = 
J r N NA — {(1, 1, 1, 1)}. The network is decomposable w.r.t. the node-arc formulation. Its 
possible decompositions are shown in Figure 



l 1 




i i 




i i 



Figure 5. Decompositions of Figure 0] with demands di = 1, i = 1,2 



In contrast to this, the network in Figure [1] is irreducible w.r.t. the path formulation. 
With the notation of Figure |5J the set of solutions of the multicommodity flow problem 
w.r.t. the path formulation is: JF P = TN F — {(1, 0, 0, 1), (0, 1, 1, 0)}. There is no non- 
trivial decomposition of the network with the set of solutions of the multicommodity flow 
problem decomposing according to it. 

Nevertheless we have a connection between irreducibility w.r.t. the node-arc 
formulation and irreducibility w.r.t. the path formulation. This connection will be 
formulated and proved in Thcorem l5.14l 
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5.5. Irreducible networks and atomic fibers in the path formulation. Now 

we want to explore the connection between irreducible networks and atomic fibers. 
We recall that a network is referred to be decomposable w.r.t. the path formulation 
for a given vector d of demand and a given vector c of capacity, if there is a 
decomposition g?i,g?2 and ci,C2 of these vectors, such that the feasible non-cyclic 
solutions of the multicommodity flow problem decompose according to it. The 
network is irreducible w.r.t. the path formulation otherwise. In fact we would like 
to prove that the irreducible networks form a certain subset of the atomic fibers, 
i.e., we want to show that they are the truncated fibers w.r.t. flow-circulations. 
We look at the multicommodity flow problem as defined in formulation H20|l . Let 
c G z|^' be a given capacity vector and d G Z^ a given demand vector. 



Let 



mm E E K U E y l P ) 

a£A 1=1 p6lli 

s.t. D-(y,s) T = (d,c) T 

y l p eZ+ VpeTL, VZ = 1,. 

s a eZ+ Va e A 



where b = (d,c) T is the right-hand side vector of 1)20(1 . denote the set of feasible 
solutions of the multicommodity flow problem. Obviously, Po.b is a fiber and the 
elements (y, s) G Po,b are in bijection with the feasible integer routings of the 
multicommodity flow problem on N = (V,A, d, c). This means: If y G T v is a 
feasible integer solution of the multicommodity flow problem on TV = (V, A, d, c) 
then (y,s) G -Po.b, where 

k 



EE 



i=i pen, 

We have the other direction, too. If (y, s) G -Pd,&, then y is a feasible integer routing 
on N = (V,A,d,c), i.e., y G T F . 

Having this bijection between the set of feasible integer solutions of the multi- 
commodity flow problem !F P and the elements in the fiber Po.b we will now establish 
the connection between the non-cyclic solutions of the multicommodity flow prob- 
lem and truncated fibers. Let L = {71, . . .7,,} denote the generating set of flow 
circulations on digraph G = (V,A) w.r.t. the path formulation. Now, we consider 
tr r (PD >6 ). 

Lemma 5.10. The elements of the truncated fiber trr(Po, b) ar ^ i n bijection with 
the non-cyclic solutions of the multicommodity flow problem J-N F . 

Proof. Suppose y G J-N F . Then (y, s) G trr(PD,&) with s defined as above. The 
other direction is clear as well. □ 

Corollary 5.11. If N — (V,A,d,c) is an irreducible network w.r.t. the path for- 
mulation, then Po,b with b = (c, d) is an atomic fiber. 

Proof. This is a direct consequence of Lemma 12.101 and Lemma 15.101 □ 

5.6. Irreducible networks and atomic fibers in the node-arc formulation. 

Up to now we have looked at the atomic fibers of the multicommodity flow problem 
with respect to the path formulation. In this subsection we will treat the atomic 
fibers of the multicommodity flow problem with respect to the node-arc formulation. 
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Consider the the multicommodity flow problem w.r.t. the node-arc formulation 
as defined in ifM)) : 

min K T f 

s.t. C-{f,s) T = (c,d,0) T d 
fls a eT+ VaeA, W = l,...,fc 

We look at the fibers of the following form: 

P c , b := {(/, s) e Z*" |A| x l}^ : C • (/, s) T = 6 T }. 

If TV = (V, A, d, c) is a network, then there is a corresponding non-empty fiber Pc,b 
with b = (c, d, 0). As for the path formulation the elements in the fiber correspond 
to the solutions of the multicommodity flow problem on N and vice-versa. 

But contrarily to the path formulation not all fibers of the matrix C correspond 
to networks on the digraph G = (V, A). Only a subset of the fibers corresponds 
to networks on the digraph G. This subset consists of those fibers with right hand 
sides b of the form b = (c, d, 0) where c e z!^', d S Z+ and the zero-vector. On 
the other hand all fibers with right hand sides of this form correspond to networks 
of the digraph G. This means that we are dealing with a monoid M of right-hand 
sides b, which can be defined as 

M := {(c,d,0) : c G z!^' , d 6 Z^_, G Z fe - (|v| - 2) }. (22) 

As in the previous subsection, we want to install a connection between the 
fibers which are atomic w.r.t. M and the irreducible networks. Therefore, let 
A = {5\, . . . ,S m } be the set of generators of flow-circulations on the digraph 
G = (V,A) w.r.t. the node-arc formulation. Then, we are able to install the con- 
nection between irreducible networks and the fibers which are atomic w.r.t. M . 

Lemma 5.12. Let G = (V,A) be a digraph, c 6 Z+ 1 ' the capacity and d E l\ the 
demand. N = (V, A, d, c) is an irreducible network if and only if the corresponding 
truncated fiber tr^(Pc^) with b = (c, d, 0) is indecomposable w.r.t. the monoid M 
in l|2l)l. 

Proof. The proof is analogous to the proof of Lemma T5. 101 □ 

Corollary 5.13. If N = (V,A, d, c) is an irreducible network w.r.t. the node-arc 
formulation, the Pc,b is an atomic fiber w.r.t. M, where b = (c,d,0). 

Proof. This is a consequence of Lemma 12.101 and of Lemma 15.121 □ 

It remains to show the connection between irreducible networks in the node-arc 
and irreducible networks in the path formulation. 

Theorem 5.14. A network which is irreducible w.r.t. the node-arc formulation is 
irreducible w.r.t. the path formulation. 

Proof. The assertion of the theorem is equivalent to the following assertion: A net- 
work which is decomposable w.r.t. the path formulation is decomposable w.r.t. the 
node- arc formulation. But this assertion is a consequence of Lemma 13.91 We have 
seen in section T5.4I that the node-arc formulation is a projection of the path-cycle 
formulation. If we denote A = {Ai,...,At} the generators of flow-circulations 
w.r.t. the path-cycle formulation, then we have: T C A and we obtain: 

A network which is irreducible w.r.t. the node-arc formulation is irreducible 
w.r.t. the path-cycle formulation. But the non-cyclic solutions of the multicom- 
modity flow problem w.r.t. the path-cycle formulation are in bijection with the 
non-cyclic solutions of the multicommodity flow problem w.r.t. the path formula- 
tion. In particular: The irreducible networks in both formulations are in bijection 
and the assertion of the theorem follows. □ 
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Figure 6. The digraph G = (V, A) for one commodity 



5.7. Formulating the network design problem as an integer linear pro- 
gram: An example. We have seen in the previous subsections that the network 
design problem for integer flows fits in our framework of the integer Minkowski 
programs because there we want to ensure that the set of solutions of the multi- 
commodity flow problem for the optimal capacity c is not empty. We have also seen 
that there is a finite generating system for the family of feasible solutions sets of the 
multicommodity flow problem. Even though there arc no Minkowski-additive and 
antitonc functionals involved in the integer Minkowski programs corresponding to 
the network design problem wc may reformulate it as an integer linear program. As 
this reformulation introduces an integer variable for each irreducible network the re- 
formulation process generates an "extended version" of the path formulation. Note 
that a vector b = (c, d, 0) T , where d = ej, i £ {1, . . . , fc}, and c = x the incidence- 
vector of a path routing d — a is necessarily contained in the set of irreducible 
networks. Therefore, all path variables y l p are contained in the reformulation as an 
integer linear program. 

We will illustrate the reformulation process by an example. 

Example 5.15. Consider the digraph G = (V,A) with three arcs and three nodes in 
Figure |S| For given d £ Z+ and given c £ Z+ the node-arc formulation for the digraph in 
Figure El looks as follows: 

fi+ si = Ci 

h + s 2 = c 2 

/s+ s 3 = c 3 (23) 

fi + h = d 

-h+h =0 

Let C denote the matrix corresponding to the left-hand side of the multicommodity flow 
problem 12311 . Then the condition of finding a feasible integer flow translates to the 
following condition: 

Pc, b = {(/, s) € Z 6 + : C ■ (/, sf = b} # 0, 

where 6 = (c, d, 0) T . Now the network design problem for the digraph G = (V, A) for given 
demand d £ Z + and vector of costs w € K+ may be formulated as an integer Minkowski 
program: 

T 

mm w c 
s.t. d — d 

(24) 

c ez+,de z+ 

Figure shows the fibers of the matrix C which are atomic w.r.t. the monoid M = 
{(c, d, 0) G Z+j. As there are no circulations possible on the digraph G, these atomic 
fibers correspond to the irreducible networks with underlying digraph G. 

Now we will reformulate the network design problem with respect to the digraph G 
and prescribed demand d £ Z+ as an integer linear program. Therefore, we introduce 
a variable A; £ Z+ for all irreducible networks with underlying digraph G. Then the 
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Figure 7. The irreducible networks with underlying digraph G = (V, A) 

reformulated network design problem is the following one: 

min wi Ai + W2 An + W3 Am + Wi Aiv + (1112 + wg)Xv + (wi + W2 + ws)\\n 
s.t. Aiv + Av + Avi = d 

Ai, . . . , Avi £ Z+. 

6. Design of survivable networks 

An optimal solution of the basic network design problem might give a capacity 
installation that admits only one feasible routing. Clearly, if a link or a switching 
node fails, not all demand can be routed through the remaining network. In this 
section, we will consider various notions of survivability of networks. The two basic 
notions are the arc survivability and the node survivability. Ensuring arc surviv- 
ability means to compute a capacity installation such that there exists a feasible 
routing of the demand even if an arbitrary link of the network fails. If a node v 
of the network fails, then all arcs incident with node v fail simultaneously and the 
commodities incident with node v, i.e., having sink or source equal to v, fail simulta- 
neously. Ensuring node survivability means to compute a capacity installation such 
that there exists a feasible routing of the remaining demand even if an arbitrary 
node of the network fails. 

The two notions can be studied in a unified and generalized failure model. We 
refer to the failure model in IPioro and Medhil <l2004) : 

Definition 6.1. A failure model is defined as a finite set S of failure states of the 
network. Each failure state a G X is characterized by a vector of arc-availability 
coefficients a a — (ax a , . . . , auu), < on t „ < 1, and a vector of demand coefficients 
Xa = (Xio- • • • i Xka); < Xjcr < !■ The arc- availability coefficients on, a represent 
the proportion of the capacity of arc i which is available in failure state a G E. The 
demand coefficients Xja denote the demand of commodity j which must be satisfied 
in failure state a G S. 

It is obvious that an arc-availability coefficient on M is equal to corresponds to 
a total failure of arc i, whereas < aj l<T < 1 corresponds to a partial failure of this 
arc and cti_ a = 1 means that there is no failure of arc i. 

In this section we will consider several particular cases for the arc-availability 
and the demand coefficients: 
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(1) x CT = landa CT G{0,l}l A l 

This case models a scenario where a single arc or several arcs a £ A fail and 
the original demand di must be satisfied for all commodities I £ {1, ... ,k}. 

(2) X a S {0, l} fc and ot„ £ {0, 1} |A| with X t* = if and only if 3v £ V„ : v £ 
{si, U} and a aa = if and only if a £ S(v) for v £V a , 

where V a denotes the set of nodes failing in failure state a. 

(3) Xo- = 1) < a.i j(J < 1 and a.j a £ {0, 1} for j ^ i. 

This situation corresponds to a partial failure of arc i and either no failure 
or total failure of arcs j ^ i. 

With the help of these "pure" scenarios one may also model mixes of them, e.g., 
simultaneous failures of several nodes, total failures of several arcs and a partial 
failure of a single arc. Note that within our consideration of arc survivability of 
networks we always mean arc survivability w.r.t. complete rerouting. This means 
that there are no routing restrictions imposed in the failure case: The no-fault and 
the fault routing may be completely unrelated. 

6.1. Arc survivability w.r.t. complete rerouting. In this section we will con- 
centrate on the survivability of networks w.r.t. total or partial failures of networks. 
The property of a network to be arc-survivable in case of failure of an arc will be 
modeled via Minkowski-additive and antitone functionals depending on the the set 
of feasible solutions of the multicommodity flow problem for fixed capacity, demand 
and a particular arc of the network. These functionals will be called arc- survivability 
functionals. Consider the functionals below: here we define them for the node-arc 
formulation of the network design problem. But it is possible to define them in an 
analogous way for the path formulation. In the following, the vector b will always 
denote b = (c, d, 0). 

A network TV = (V, A, d, c) is survivable w.r.t. complete rerouting in case of total 
failure of arc a £ A if and only if g a (Pc,b) = 0. 

Lemma 6.2. Let a £ A be fixed. The functional g a (Pc,b) is Minkowski-additive, 
i.e., ifN = (V,A,d,c) is decomposable into N 1 = (V, A, d 1 ,^) and N 2 = (y,A,d 2 ,c 2 ), 
then we have g a {Pc,b) = ga{Pc.b 1 )+ya(Pc,b 2 ) ■ Furthermore, g a {Pc.b) is an antitone 
and non-negative functional. 

Proof. Let TV = (V, A, d, c) be a network and a £ A be an arc. Let /* be the 
solution of the multicommodity flow problem H20JI on N with minimum aggregated 
flow on arc a. Then we have: 

k 

9a(Pc,b)=Y,(na- 
1=1 

W.l.o.g. we may assume, that /* £ TN NA , i.e., /* admits no circulation of flow. 

Now let N = Ni + N 2 , implying that f*=fi + f 2 with f t £ FNf A , i = 1,2. 
This yields immediately that the solutions of the multicommodity flow problem 
(I20|l on Ni, fi, i — 1,2, admit each minimum aggregated flow on arc a. Therefore 
we have: 

g a (Pc,b) = 9a(Pc,b l ) +9a(Pc,b 2 )- 
This implies finally the Minkowski-additivity of the arc-survivability functionals 
g a (Pc,b)- 

The arc-survivability functional g a {Pc,b) is antitone as the minimum-function is 
antitone. Its non-negativity comes from the non-negativity of feasible flows. □ 



INTEGER MINKOWSKI PROGRAMS 



23 



Since g a {Pc,b) is Minkowski-additive and antitone we can apply Theorem 12.61 
We set 8 := and add the inequalities 

t 

i=l 

to the reformulation of the network design problem as integer linear program to 
ensure survivability w.r.t. complete rerouting in case of failure of an arc. As the 
functionals g a (Pc,b) are non-negative we can replace these inequalities by the fol- 
lowing equalities: 

t 

9a{Pc,b i )=° VaeA 

i=l 

and arrive at the following reformulation of the network design problem under arc 
survivability constraints as integer linear program: 

min W T {jz\ Cl ) 
i=i 

t 

s.t. J2 Xidi = d 

«f (26) 
£ 9a(Pc,bi) = Vaei 

tci,di,0) T = h GF(C) Vi = l,...,t 

A, G Z + Vi = l,...,t 

Remark 6.3. We remark that this integer linear program is infeasible if the topol- 
ogy of the network does not support a survivable installation of capacities. This case 
will occur if there is a commodity pair that does not have two arc-disjoint paths. 

We will illustrate this reformulation process with the help of an example. 

Example 6.4. We resume Example 15 . 1 5 1 and incorporate the arc-survivability functionals 
in our reformulation. Consider Figure |H1 The values of the arc-survivability functionals 
on the arcs of the irreducible networks are marked in bold. 




Figure 8. The values of the arc-survivability functionals 
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Now we may consider the reformulation of the network design problem under arc- 
survivability constraints. For our example this reformulation is the following one: 

min wi Xi+ W2 Xn+ «)3 Ami Wi Aiv+ (11)2 + W3) Xy+ (wi + W2 + W3) Avi 
s.t. Aiv+ Av+ Avi = d 

0Ai+ OA11+ 0A m + 1Aiv+ 0A V + 0A V i = 

0Ai+ OA11+ 0A m + 0Aiv+ 1A V + 0A V i = 

0Ai+ OA11+ 0A m + 0Aiv+ 1A V + 0A V i = 

Ai, . . . , Avi G Z+. 

Remark 6.5. Note that if a network is arc- surviv able w.r.t. the failure of an arc 
it will decompose into survivable, irreducible summands. This implies that all vari- 
ables Xi for non- surviv able, irreducible networks JVj are fixed to in the reformula- 
tion as integer linear program of the survivable network design problem. 

6.2. Partial failures. Up to now, we have considered total failures of arcs. But it 
is possible to model so-called par tial failures of arc s . We want to clarify this with 
the help of the failure model of IPioro and Medhil l|2004|) . The above discussion 
has enabled us to model the |S| = \A\ failure states with x<r a = 1, aw a = and 
a eaa = 1 for all a =/= e G A. The above model guarantees survivability of the 
network in each of these failure states. But we may model partial failures of arcs 
in an analogous way. To this aim let E with |E| = \A\ be the index set of the 
failure states with Xa a = 1, < a aa - a < 1 and a e( ? a = 1 for all a ^= e G A. Then 
incorporating the following inequalities to our model guarantees survivability of the 
network in all of these failure states: 

g a (c, d) < (27) 

These constraints fit in our framework of integer Minkowski programs: in section 
12.11 we have seen that we may incorporate non- linear constraints 

k(A(z)) + h T z < 5 

with k Minkowski-additive and antitone. We already know that the arc-survivability 
functionals g a {Pc,b) are Minkowski-additive, antitone and non-negative. Now let 
h a := — a aMa e a with e a G zli 4 x )\ v \ ^} 1C unit-vector admitting a 1 in the a-th 
position and let r := 0. Then the condition 

ga(Pc,b) + h T b<T 

is equivalent to inequality (|27|l . 

6.3. Simultaneous failures. Up to now, we have considered the case of surviv- 
ability w.r.t. complete rerouting in case of failure of only one arc. Our approach to 
model survivability of a network in case of failure of an arc can easily be modified to 
model survivability w.r.t. complete rerouting in ca se of simultaneous failure of up 
to q arcs. In the model of lPioro and Medhl l|200-j) this means: Let £ be an index 
set of failure states with |S| < 00. \ a = 1 for all a G S and a a G {0, l}'" 4 '. Now 
let a & S be fixed and let a±, ■ . . , a q be the set of arcs with a 0i;(T = 0. Consider the 
Minkowski-additive and antitone functionals presented below: 

g ai ,..., aq (Pc, b ) ■■= b f E (fl +fL + --- + fa q )) ■ 

We will call functions of the above type arc-survivability functionals for simulta- 
neous failure of arcs a\, . . . , a q . The Minkowski- additivity of these functionals can 
be shown analogously to the proof of Lemma f6. 21 Antitonicity and non- negativity 
of the function follow as well. As before, a network TV = (V, A, c, d) is survivable 
w.r.t. complete rerouting in case of simultaneous failure of arcs ax,...,a,q if and 
only if g ai ,...,a q {Pc,b) = 0. Applying Theorem 12.61 and adding these equalities to 
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the formulation of the network design problem as integer Minkowski program for all 
failure states a £ E ensures arc survivability of the network. As before survivability 
can be guaranteed by adding linear equalities to the formulation as integer linear 
program. 

6.4. Node survivability w.r.t. complete rerouting. The arc-survivability func- 
tionals for simultaneous failure of q arcs will allow us to model node survivability 
of a network. If a node v £ V fails then all arcs incident with this node fail si- 
multaneously. Additionally the demand of all commodities incident with v fail as 
well. This implies in particular that the load of the network is lower in the failure 
case! Therefore we have two different cases: Either node v is the source or sink of a 
commodity or it is not. We consider a transformation of our digraph (sec Figure|5J) 
which enables us to model node survivability of a network via arc-survivability 
functionals. For each commodity I we introduce an auxiliary arc ai = (si,ti). Our 
digraph now consists of the set of nodes V and the disjoint union of ordinary arcs 
A and auxiliary arcs A, i.e., G = (V, A OA). Each arc a; £ A has a capacity of 
c a , = di. If a node v fails we want the commodities incident with this node to be 
routed via the auxiliary arcs a £ S(v). 




FIGURE 9. Transforming the digraph G = (V, A) to G' = (V, A U A) 

Now we will model node survivability of a network. To this aim let v £ V be 
fixed. We wa nt to model survivability of the network in case of failure of node v. In 
the model of IPioro and Medhl l|2004|) the corresponding failure state would be the 
following one: x<j £ {0, l} fc with Xal = if and only if v £ {si,U} and a. a £ {0, 
with a aa = if and only if a £ 5(v). Now we want to model survivability of the 
network in this failure state. First of all we have to ensure the existence of a no-fault 
routing which does not use any of the auxiliary arcs and satisfies all the demands. 
This means that the minimum simultaneous flow on all auxiliary arcs is equal to 0: 

<7« ,l'r>-< "• (28) 

Additionally one has to ensure that there exists a feasible vector of flow in case 
of failure of a node v. Then the following constraint ensures survivability of the 
network in case of failure of arc v. 

(Pc,b) = (29) 

aeA\8(v) 

With the help of these auxiliary arcs and the arc-survivability functional for si- 
multaneous failure of a set of arcs it is possible to model node survivability of a 
network by adding the above constraints to our formulation of the network design 
problem as integer Minkowski program. As shown in the previous section it is possi- 
ble to guarantee node survivability by adding linear constraints to the reformulation 
as integer linear program. 

It is clear that one may model more complex failure situations by combining and 
adapting the above scenarios. 
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Conclusions 

In a forthcoming paper l|Eisenschmidt. Konne. and Laugierll2006btL we will present 
algorithms to compute atomic fibers for network design problems and computa- 
tional results. We remark that we expect that the integer linear programs arising 
from the reformulation method will be extremely large. It is therefore desirable 
to devise a method for computing atomic fibers on the fly, and to use it within 
a specialized branch-cut-and-price algorithm. However, algorithmic questions like 
this are beyond the scope of this paper; they will be the topic of future research. 
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